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ABSTRACT 

Buoy  data  provided  clear  evidence  of  mixed  layer  deepening  and  an 
internal  wave  caused  by  Hurricane  Eloise , Septeober  1975 . Logarlt±nilc 
tenperature  profiles  below  an  isotheimal  mixed  layer  were  assumed  and 
used  to  model  thexmocUne  oscillation  and  heat  budget  calculation  as 
influenced  by  Eloise  over  a 21-day  period.  Resiilts  show  tbat  prior  to 
the  arrival  of  Eloise  at  the  buoy,  the  average  mixed  layer  depth  was 
about  33m.  As  the  winds  increased  due  to  hiirrlcane  approach,  the  mixed 
layer  deepened  steadily  to  about  42m  before  upwelling  to  approadaately 
22m.  The  thermocline  then  underwent  three  distinctly  large  os<d.llatlons 
of  inertial  periodicity,  vhile  the  mixed  layer  ccxitinued  to  deepen.  The 
post-storm  average  mixed  layer  depth  was  abcxit  S2m.  Val\jes  of  mixed 
layer  depth  were  concluded  to  be  accurate  to  within  2&t.  Vertical  velo- 
cities, calculated  first  by  assuning  zero  horizontal  taiperature  advec- 
tlon  in  die  material  derivative  ecjuatlon  and  second  by  finding  the  mass 
transport  necessary  to  balance  the  heat  budget,  shew  that  in  the  upper 
500m  of  the  water  colum  downward  vertical  motion  of  Im/hr  or  less  pre- 
vailed during  storm  approadi,  followed  by  upward  vertical  velcxdty  as 
great  as  5.35m/hr  during  the  12  hr  iamediately  following  hurricane  pas- 
sage followed  by  downrard  vertical  velcxdty  during  the  large  thermocline 
osdllatlons . 
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I.  INIRODUCnON 


By  a stroke  of  Ivick,  the  eye  of  Hurricane  Eloise  on  23  Septenber  1975 
passed  over  the  National  Oceanic  and  Atmospheric  Administration's  EB-10 — 
a forty- foot  diameter  data  buoy  andiored  in  the  Gulf  of  Mexico.  Numerous 
Instnxoents,  both  on  the  buoy  Itself  and  attached  to  the  line  anchoring 
the  buoy  to  the  floor  of  the  Gulf,  gathered  atmospheric  and  oceanographic 
data  througjhout  the  three  days  before  and  18  days  after  hurricane  pas- 
sage. The  purpose  of  this  thesis  is  to  study  the  response  of  the  upper 
ocean  to  Hurricane  Eloise  and  to  report  the  resiolts  obtained. 

The  medianisms  diat  produce  changes  in  ocean  thermal  structure  in  the 
wake  of  severe  tropical  cyclones  have  been  studied  for  seme  time.  In 
addition  to  the  heat  loss  to  the  atmosphere , the  oceanic  processes  of 
vertical  and  horizontal  advection  plus  turbulent  mixing  at  the  top  of  the 
thermocline  have  beai  shewn  to  contribute  to  cooling  of  the  upper  layers 
in  the  ocean.  Jordan  [1964]  reasoned  that  for  typical  mixed  layer  depths, 
large  tesperature  decreases  could  not  be  due  to  the  heat  loss  to  the 
atmosphere  but  mist  originate  from  the  thermocline  layers  througji  verti- 
cal fluxes.  Leipper  [1967]  reported  ^diat  is  probably  the  best-knewn 
survey  of  oceanic  thermal  response.  Black  and  Mallinger  [1972]  coobined 
airborne  expendable  bathythermograph  and  conventional  data  to  study  the 
effects  of  hurricane  Ginger.  Althou^  a comnon  shortcoming  of  the  obser- 
vational studies  has  been  the  sparsity  of  coincident  before-after 
soundings  allowing  determination  of  the  changes  due  to  the  hurricane, 
Fedorov  [1973]  used  before  and  after  soundings  at  ocean  weedier  station 
(OUS)  Tango  to  calculate  the  teoperature  changes  due  to  the  passage  of 


14  tTphoons,  Elsberry  et  al.  [1974]  developed  an  enpirical  hurricane 
boundary  layer  model  for  hurricane-ocean  interaction  studies.  This 
model,  vAiidi  was  used  recently  to  drive  a mixed  layer  model  of  the  upper 
ocean  [Elsberry  et  al. , 1976]  to  sumlate  the  dieimal  response  induced 
by  a hurricane,  enphasized  the  role  of  vertical  fluxes  by  incorporating 
a mixing- layer  model  similar  to  Kraus  and  Turner  [1967]  and  Dennan  [1973] . 
In  tdiis  thesis,  logarithmic  temperatinre  profiles  as  suggested  by  Tully 
[1953]  were  used  to  estimate  the  temperature  structure.  These  profiles 
were  used  to  predict  the  depth  of  the  mixed  layer  and  to  balance  the 
heat  bxidget  in  the  wake  of  the  storm.  Oscillations  of  the  thermocline 
predicted  by  Geisler  [1970],  docunaented  by  Black  [see  Sheets  1974]  and 
modeled  by  Grigsby  [1975]  were  observed.  The  heat  budget  was  calculated 
producing  the  horizontal  and  vertical  velocities  necessary  to  account 
for  t±ie  observed  temperature  changes. 

II.  PROCEDURE 

A.  DATA 

The  buoy  EB-10,  located  at  27“28'N,  88“01’W,  in  D13  fathons  of 
water,  gathered  atmospheric  data  plus  the  temperature  and  hydrostatic 
pressure  of  the  water  near  the  surface  and  at  three  sihsurface  depths 
in  the  Gulf  of  Mexico  during  the  passage  of  Hurricane  Eloise.  These 
data  were  published  by  the  lOAA  Data  Buoy  Office  in  a report  titled 
Data  Report;  Buoy  Observations  During  Hurricane  Eloise  (Septenher  19 
to  October  11,  1975,  Novenfcer  7,  1975)  (fei thee  and  Johnson,  1975). 

B.  TEMPERATURE  PROFILES 

To  study  die  nature  and  extent  of  influence  of  a hurricane  on  the 
ocean,  having  the  teoperature  of  the  sea  water  before,  during  and  after 
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hurricane  passage  at  all  depths  would  be  optimal.  But  the  system  avail- 
able provided  temperatures  with  a maxinimi  frequency  of  only  once  per  hour 
and  recorded  them  only  near  the  surface  and  at  three  sub-surface  depths — 
approodmately  50,  200  and  500  meters.  These  are  only  nominal  depths  be- 
cause of  the  vertical  wandering  of  10-12  meters  due  to  periodic  pulling 
on  the  mooring  line  by  the  drifting  buoy  in  the  hi^  winds  and  accompa- 
nying waves . The  near-surface  temperature  value  was  measured  at  a depth 
of  two  meters  and  will  be  referred  to  throughout  this  thesis  as  the 
surface  temperature. 

The  first  goal  of  this  thesis  was  to  estimate  the  temperature  at  all 
depths  for  any  observational  time  kncwing  the  teiqperature  at  only  the 
surface  and  three  subsurface  depths . Estimates  of  the  above  type  were 
maH«»  for  170  different  observation  times  throughout  the  21-day  period 
bracketing  the  passage  of  the  storm  over  EB-10 . A temperature  structure 
of  the  water  was  also  estimated  by  the  NQAA  Data  Buoy  Office  0®B0) 
[Withee  and  Johnson,  1975,  page  B-30] . A portion  of  NDBO's  plot  is  shewn 
in  Figure  1;  the  surface  and  three  subsurface  temperature  values  are 
written  for  eadi  observation  time  at  the  depth  of  the  measuirements . The 
isotherms  shown  in  Figure  1 were  drawn  assuming  a linear  tesperature 
change  with  depth.  Blade  and  Withee  [1976]  and  Price  [1977]  used  these 
same  data,  and  estimated  the  temperature  change  with  depth  for  all  times 
with  gradients  based  on  a single  AXET  dropped  near  EB-10  on  the  day  prior 
to  Eloise's  passage  over  the  buoy.  In  this  thesis  the  positions  of  iso- 
theniB,  as  well  as  the  depth  of  the  theimocline , were  calculated  at  each 
of  the  170  observation  Hnip«  assuming  a logarithmic  temperature  change 
with  increasing  pressure  defined  by  the  three  observed  subsurface  tenper- 
ature  valujes.  Eadi  logarithmic  curve  was  then  matdued  to  a vertical  line 
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representing  an  Isothermal  mixed  layer  having  teoperature  equal  to  the 
surface  temperature  observed  at  the  respective  observation  time.  This 
procedure  assunes  the  three  svjhsurface  temperatures  were  within  the 
thermocllne.  After  passage  of  Eloise,  the  nominal  50  m observations 
were  occasionally  equal  to  the  surface  temperature,  Indicating  that  the 
mixed  layer  was  greater  than  50  m.  In  these  cases  only  the  two  lowest 
s\jbs\jrface  observations  are  available  to  define  the  logarithmic  taopera- 
ture  profile. 


1.  Procedure  For  Calculating  The  Logarithmic  Temperature  Profile 
Below  The  Mixed  Layer  At  Any  Givm  Observation  Time 


At  0900  GMT,  21  Septedber  1975,  the  data  provided  by  EB-10  indi 


cated  that  the  three  following  subsurface  temperatures  existed  at  the 
given  pressxxres  : 


^ob 

P 

In  p 

25.784“C 

6. 42 7kg/ cm^ 

1.8605 

15.677 

22.986 

3.1349 

9.055 

52.148 

3.9541 

fdiere  is  the  observed  tasperature  and  p is  pressure.  A least- 

squares  linear  regression  was  calculated  using  the  natural  logarldm  of 
pressure  (In  p)  versus  ^ob  • 

Figure  2 shows  graphically  the  best-fit  curve  throu^  the  three 
sets  of  values  being  considered.  The  slope  and  intercept  of  the  regnis- 
sicn  curve  for  this  case  were  calculated  as  -0.125  and  5.092  respec- 
tively. That  is,  Che  equation  celating  pressure  and  cemperatuce  is 

In  p - -0.125T  + 5.092  . (X) 

A similar  equation  was  calculated  for  eadi  of  the  170  observa- 
tion times  using  the  Naval  Postgraduate  School  IBH  360  computer 
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FIGURE  2.  Example  best-fit  linear  regression  curve  obtained 
from  sxi>surface  tenperature  and  pressure  values.  This  example 
Is  for  die  observation  time  of  0900  GMT,  21  September  1975. 
Curve  has  a slope  of  -0.125  and  an  Intercept  of  5.092. 


subroutliie  to  calculate  linear  regressions.  The  slopes  m and  inter- 
cepts b specify  the  thern)cx:line  ten|)erature  profiles  dtrou^iout  the 
21  days  \nder  consideration.  The  general  eqviation  is 

In  p - off  + b . Cla) 


To  estimate  the  accuracy  of  equation  G)  > Che  three  values  of 
the  logarittDi  of  pressure  used  in  the  regression  can  be  used  in  equa- 
tion G)  to  calculate  the  regression  temperatures  GjJ  • The  following 
infozmation  was  foizid: 


Error  - 


2 

p Ocg/cm  ) 

Tob  C"C) 

Tj.C“C) 

Tr  - C-C) 

6.427 

25.819 

25.804 

0.015 

22.986 

15.677 

15.627 

-0.050 

52.148 

9.055 

9.085 

0.030 

The  largest  error  in  temperature  for  this  observation  time  is  -0.05*C, 
and  is  associated  with  the  middle  of  ^le  three  stjbsurface  temperature 
measurements— a diaracteristic  of  drawing  the  best  linear  fit  throu^ 
the  three  points.  According  to  WBO,  the  accuracy  of  the  temperature 
sensors  was  O.OS^C  for  the  stirface  and  two  deepest  subsurface  sensors, 
and  0.2*C  for  the  sensor  located  near  the  50  m depth.  Consequently  the 
logarithmic  profile  tedmique  fit  the  data  at  the  three  levels  within 
the  expected  accuracy  of  die  temperature  observations.  A more  complete 
description  of  teiperature  profile  errors  will  be  given  later. 

2.  Procytoe  For  Osculating  The  TemperaG^  and  Depth  of  The 
Mixed  Layer  At  Any  tiiven  Cbservation  Thne 

The  next  step  in  die  procedure  was  to  matdi  the  observed  sur- 
face temperature  with  the  theimocline  profile  given  by  equation  (1) . 

The  surface  temperature  at  0900  GMT,  21  September  1975  was  observed  to 
be  28.819*C.  Ehterlng  this  value  for  temperature  in  equation  G)  gives 
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p ~ 4.A06.  Using  this  calculated  pressure  value  and  assxnlng  an  Isotiier^ 
mal  adjced  layer,  the  ccnplete  tenperature  profile  can  nov  be  drawn  fixm 
the  sxjrface  to  the  point  of  the  deepest  observation.  This  conoplete  pro- 
file for  0900  GMT,  21  Septeoher  1975,  Is  shown  in  Figure  3^.  Figure  3B 
shows  the  sane  profile,  except  that  pressure  was  plotted  linearly  vice 
logarithmically . 

Figure  3C  shows  die  profile  assxming  a linear  taiperacure  change 
in  depth  as  used  by  the  NQAA  Data  Buoy  Office  when  plotting  Figure  1. 

The  profile  in  Figure  3B  appears  to  be  a nudx  more  realistic  representa- 
tion of  conditions  observed  dally  in  the  subtropical  oceans.  The  point 
of  intersection  of  the  vertical  isotheimal  line  with  the  regression  curve 
in  Figure  3A  or  3B  Is  defined  in  this  thesis  as  the  hydrostatic  pressure 

at  the  bottom  of  the  mixed  layer,  l.e.  at  the  top  of  the  thezmocline . 

2 

Throu^out  this  thesis,  pressure  p Ocg/on  ) was  converted  to 
depth  z (m)  according  to  an  approximation  algorithm  suggested  by  M)B0 
for  the  sxhsurface  Ooore  than  AOm)  data 

z - C9.75p  - 10)  . (2) 

Hence  the  mixed  layer  depth  can  be  found  for  0900  (IfF,  21  Septaher  1975 
by  substituting  p * 4.A06  into  equation  C2)  giving  MLD  ■ 33. Qm. 

Temperature  profiles  extending  from  the  siirface  to  the  deepest 
observation,  as  well  as  the  respective  mixed-layer  depths,  were  calcu- 
lated for  all  170  observation  tines  as  in  the  example  above.  A plot  of 
MLD  vs  time  was  made  for  studying  the  motion  of  the  thezmocline  and  the 
deepening  of  the  mixed  layer  during  the  three  days  before  and  18  days 
after  the  passage  of  Hurricane  Elolse  over  EB-10.  This  plot  will  be 
presented  and  discussed  in  a later  section. 


16 


FIGURE  3.  Derived  tgqperature  profile  siirface  Co  point  of  deepest 
observadcn  for  the  obaervatlon  tine  of  0900  OfT,  21  Septeober  1975. 
Pressxsre  is  plotted  logarltimlcally  in  A and  linearly  in  B ; C shows  the 
profile  for  this  observation  tinie  assuned  by  NCAA  Data  Buoy  Office.  An 
isothextnal  oxlxtd  layer  was  assuned,  i.e.  a vertical  line  was  draNn  &cm 
sxjrface  temperature  to  point  of  incersecdcn  with  regression  curve. 


3.  Procedure  For  Calculating  The  Position  of  I80^h«»rTT» 


The  depth  z of,  e.g. , the  25”C  isotherm,  can  be  found  for 
0900  GMT,  21  Septenber  1975  by  letting  T - 25  In  equation  Q)  and 
solving  for  p and  then  substituting  for  p into  equation  (2}  • Doing 
this  gives  a depth  of  59.3m  for  the  25°  isotherm  at  this  time.  For  this 
thesis  depths  were  thus  calculated  for  the  25,  23,  21,  18,  15,  12  and  9° 
isotherms  at  all  170  observation  times . A plot  of  isotherm  depth  versus 
time  was  made  for  studying  isotherm  motion  duzring  the  three  days  before 
and  18  d^s  after  the  passage  of  Hurricane  Eloise  over  EE- 10.  This  plot 
will  be  coqpared  to  a similar  plot  made  by  NCAA's  Data  Buoy  Office  shown 
in  Figure  1 above. 

C.  CALCULATINC;  THE  HEAT  BUDGET 

The  procedure  used  in  this  thesis  to  calculate  the  heat  budget  was  a 
sinplified  one.  The  problem  was  unique  in  that  the  heat  budget  was  exa- 
mined over  a period  of  time  with  large  amplitude  oscillations  on  the 
thermocline.  The  time  change  of  heat  content  (cal/cm  ) of  a colum  of 
water  imder  the  buoy  500  meters  deep  and  one  square  centimeter  in  cross 
section  is  related  to  the  time  change  of  average  tenperature  obtained 
from  the  profiles  derived  above.  That  is, 

- eCp  ' “av  O) 

2 

where  H * heat  in  cal/on 

3 

p "■  density  of  water  ^ 1 gn/aa 
- heat  capacity  of  water  at  constant  pressure  * 1 cal/gaa/'*C 

z ••  the  height  Cob)  of  the  colmn  of  water  one  square  centimeter 
in  cross-section 

T " average  tsoperature  C*C)  of  the  coltxn  of  water,  as  calcu- 
lated using  the  logarlthnlc  regression  curves  discussed  ear- 
lier 
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In  temis  of  Cinie  derivatives  equation  (3)  becomes 


r 


i 


at 


3Tav 

at 


C3a) 


and  will  be  called  the  storage  teim.  If  the  tenperature  change  is  loea- 

sured  in  degrees  Celsius  per  hour,  the  storage  term  will  have  vnits  of 
2 

cal/cm  /hr. 

From  the  first  law  of  diemxxfynamics , two  mechanisms  were  con- 
sidered as  changing  the  heat  content  of  the  colum — the  net  flux  of 
latent,  sensible  and  radiative  heat  transfer  across  the  sea  surface,  and 
heat  flux  by  three  dimensional  advection.  The  formulae  used  [see  Husby 
and  Seckel,  1975]  for  obtaining  the  net  flux  of  latent,  sensible  and 
radiative  heat  (0  ) was 


^ - Qs  - ^ - 9e  - Qc 


viiere  , the  solar  insolation  term,  was  measured  directly  by  ^-10, 

- 1.1^0'^ (273. 1&4T^^ 
xC0.39-0.05e^*^2)Cl-0.6C^)  , 

cl 

0 - 3.767C,C0.98e  -e  )W  , and 
e u w a 

- 2.488CpCr^-T^)W 


e is  die  saturation  vapor  pressure  (nb)  of  the  aonosphere  at 
^ die  hel^t  of  ten  meters  calculated  using  die  Goff-Gratdi 
fomulatlon  of  die  Clausius-Cl^peyron  equation  with  the  dew^ 
point  teoperature  as  the  entering  argument 


e is  die  saturation  v^ior  pressure  (nb)  of  die  atmosphere  at 
die  sea  surface  calculated  similar  to  e except  using  sea 
surface  tenperature  as  die  entering  argi&nt 

C.  is  the  nondimenslonal  drag  coefficient  equal  to  (0.63  + 

“ 0.66*i7)*10-3,  [Snlth  aid  Baike,  1975] 

T is  the  teoperature  of  the  air  in  degrees  Celsius 

is  the  tenperature  of  the  water  in  degrees  Celsius 

W is  the  wind  speed  in  meters  per  second 
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«nd  C,  (the  fraction  of  sky  covered  by  clouds) , 

Is  0.80  from  1200  GMT  20  SEP  to  1800  GMT  22  SEP, 

1.00  from  2100  Off  22  SEP  to  1200  GMT  23  SEP. 

0.85  from  1500  GMT  23  SEP  to  1200  GMT  24  SEP, 

0.75  from  1500  Oil  24  SEP  to  1100  GKT  25  SEP,  and 

0.50  from  1200  Off  25  SEP  to  1200  GMT  26  SEP. 


2 

Herein  (cal/on  /hr)  is  called  the  surface  flux  term. 

The  500-meter  deep  colunn  was  then  divided  Into  five  layers , 
each  100  meters  deep,  and  vertical  and  horizontal  velocities  were  calcu- 
lated that  satisfied  the  heat  budget.  The  first  atten^jt  to  obtain  these 
velocities  was  the  commonly  used  assuiption  that  the  material  derivative 
of  temperature  was  zero  and  that  there  was  no  horizontal  temperature 
gradient.  Starting  with  the  equation 


3t 


0 


— + V 
3t  ''h 


• TT  + w ^ 


3z 


(5) 


setting  Vjj  • VT  “ 0 , and  solving  for  w gives 


w 


3T  , ^ 
3t  ' 3z  • 


(5a) 


Both  31/ 3t  and  3T/3z  were  directly  available  frcm  the  170  teopera- 
ture  profiles  obtained  as  discussed  earlier.  The  values  of  w were 
dius  cooputed  at  100  , 200  , 300  , 400  and  500  meters.  These  values  of  w 
were  related  to  the  vertical  mass  flow  rate  into  or  out  of  each  layer 
of  the  colvim  by  the  equation 


M,  - pw  (6) 

z 

vhere  is  directed  positive  upward  in  accordance  with  the  positive 
z axis. 
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From  die  ccnCinulty  equation  the  horizontal  mass  flux  into  or 
out  of  the  colum  was  calculdile  from  the  vertical  mass  fluxes.  The  net 
horizontal  mass  flow  must  ccopensate  for  the  difference  in  vertical  mass 
fluxes.  Let 


pVjj  dz 


vdhere  - the  horizontal  transport  of  mass  per  unit  width  of  vertical 
surface 

Vjj  ■ the  horizontal  velocity, 

Zg  and  * the  depths  in  meters  at  the  bottom  and  top  of  the  colum. 
For  inconpressib  le  flow  in  a constant -density  fluid, 


V • pVjj  dz  + 


J ^ 9z 


dz  - 0 


CoDDbined  with  equations  (6)  and  C7) , the  above  equation  gives 


- M - M . (8) 

Horizontal  mass  flow  is  defined  as  positive  if  directed  out  of  the 


colum. 

There  was  some  question  vhether  the  w's  obtained  using  equa- 
tion (5a)  could  be  accurate,  since  die  assvnption  of  zero  horizontal 
taiperature  advection  seemed  unlikely.  This  was  checked  by  using  the 
calculated  vertical  velocities  and  the  storage  and  surface  flux  tens 
in  the  heat  budget.  The  total  heat  flux  (THF)  across  all  faces  of  a 
water  colum  due  to  three-dimensional  advection  was  defined  as 


- Vp’i*  ■ ■ W*'' 
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where  Is  the  average  tesoperature  of  a lOOzn-deep  column  of  water 
with  cross-section  equal  to  one  square  centimeter,  and  Xj,  the 
tenperatures  at  the  botton  and  top  respectively  of  the  same  lOOm-deep 
colum,  and  M , M and  are  the  mass  transports  across  the 

bottom,  top  and  horizontal  faces,  respectively,  of  the  colum.  Substi- 


tuting equation  (8)  into  (9)  and  adding  the  storage  and  surface  flux  terms 
give  the  full  equation  used  in  analyzing  the  heat  budget. 

Rearranging  provides  the  equation 


\ ■ ■ s * ■ • V - ^ 


- V 


Cooblning  equation  (10)  with  equatioL.  (6)  gives 

w-Vp-  V ■ 


/p  (10a) 


Vtiile  M was  considered  to  be,  by  the  kinematic  boundary  con- 

dition,  equal  to  zero  at  the  sea  surface,  was  considered  to  be 
equal  to  zero  at  all  depdis  except  the  surface.  Because  the  M , M 

and  are  meoi  quantities  averaged  over  the  time  increment  between 
observetiona,  the  contribution  of  turbulent  motions  on  smaller  time 
scales  is  not  considered. 

As  an  example  consider  the  problem  of  calculating  the  vertical 
velocity  of  100  m as  given,  first,  by  equation  (5a)  and  second  by  equa- 
tion (10a)  for  the  observation  time  of  1500  GMT,  23  Scpteober  1975: 


1.  Calculation  of  Vertical  Velocities  Assuning  Zero  Horizontal 
temperature  Advectlon 

The  values  of  slope  m , intercept  b,  pressure  at  the  bottom 
of  the  mixed  layer  and  MLD  in  the  following  table  were  obtained 

from  the  respective  regression  curves  plus  equations  (la)  and  (2) , v*iile 
the  surface  tenperature  was  obtained  from  the  pihlished  btjoy  data 

for  23  Septenber  1975. 


'^ob 

m 

b 

MLD 

^sfc 

^lOOm 

car 

(kg/on^) 

(m) 

(“C) 

(TC) 

1200 

-0.152 

5.345 

3.331 

22.47 

27.226 

19.205 

1500 

-0.153 

5.334 

3.280 

21.98 

27.166 

19.072 

1800 

-0.145 

5.253 

3.733 

26.40 

27.166 

19.532 

The  mmerator  of  equation  (5a) , 

refers 

to  a single  level,  say 

100m,  and  is  evaluated  at  1500  GMT,  using  a centered  finite  difference, 
as 

(19.532  - 19.205)*C/6hr. 


The  denominator  in  equation  (5a)  was  approximated  by  using  the  1500  (211 
regression  teoperature  profile  over  an  increment  one  meter  above  and  belcw 
the  100  meter  level.  Thence  the  calculated  value  of  vertical  velocity 
w was  -0.9Im/hr.  Similar  calculations  produced  the  values  for  this 
method's  vertical  velocities  at  eadi  of  the  deeper  levels.  The  left- 
hand  set  of  arrows  in  Figure  4 shows  die  values  calculated. 

1 

2.  (>1<^tTai-ion  of  Vertical  Velocities  Not  Assuning  Zero  Horizontal  i 

Tanperature  Advectlon  " 

To  solve  for  the  mess  transport  at  the  100  meter  level  for  15(X) 

(2Cr  23  Septenber  1975  xising  equation  CTO) , we  must  calculate  the  aver- 
age tenperature  in  the  upper  100m  of  the  ocean  in  the  vicinity  of  the 
buoy  for  each  of  the  three  observation  times  listed  in  paragraph  1 
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w (m/hr) 


Wx/p  (m/hr) 


FIGURE  4.  Vertical  velocities  at  100m  depth  increments  in  the  left 
coliim  viere  obtained  using  equation  (5a) , l.e.  assiming  zero  hori- 
zontal tenperature  advectlon.  Values  of  M /p  In  rl^t  colum  were 
obtained  using  equation  (10a) , l.e.  not  ass&dng  zero  horizontal 
tenperature  advectlcn.  These  vertical  velocities  are  for  the  obser- 
vation time  of  1500  (KT.  23  September  1975. 


0 


above.  Assume  the  general  equation  is 


^a\p0^100m  “ WIXT^fc)  + C100.MLD)(T3^^ij(^1/100  (11) 


The  average  tecperature  for  any  layer  of  water  having  a logarith- 
ndc  temperature  profile  and  lying  between  press\xre  surfaces  Pg  and 
pkp  is  given  by 


^av^pj,  ^ 


1 

Pg-Pr 


t 


T(p)dp 


(12) 


Ccnbining  equation  (la)  with  (13)  gives 


In  p-b 


m 


^ "*^3  Pb'^T 

m(pj.-Pg)[pg(ln  Pg-pg-bpg)-  ( Li  p.j.-p.j,-bp^)  ] 


•)  dp 


(13) 


Equation  (13)  can  be  used  to  get  the  average  ten?>eracure  frcin 
the  bottom  of  die  mixed  layer  to  the  100m  depth,  ^Aiich  for  1200  QiT, 

23  September  turns  out  to  be  22.422'’C.  This  value  can  then  be  used  in 
equation  (11)  to  get  die  average  temperature  for  the  layer  O-lOOm,  vbich 
is  23.501*’C  in  this  example.  In  an  analogous  way,  the  average  tempera- 
ture for  die  upper  100m  of  water  at  1800  CUT,  23  September  1975  is  found 
to  be  23.8A9®C.  Using  this  information  and  equation  (3a)  finite' dif- 
ferenced about  1500  GMT,  23  September  1975,  the  change  in  heat  content 
of  the  100m  colum  is  calculated  as 


(M) 

''  At  ''l500GMr 
• lOOOOon  ■ 


" '^p^^S^1200-1800GMr 
C23.849-^..5pi>-,C.  . 580cal/cm2/hr 
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2 

The  svnface  flux  tern  calculated  from  equation  (4)  was  45  cal/om  /hr. 
The  only  remaining  variables  needed  to  solve  for  the  mass  transport  at 

100m  are  (T^^g j^gg^  and  for  1500  GMT.  These  are  available 

from,  respectively,  equation  (11)  and  the  regression  curve  confuted  for 
1500  GMT  and  are  equal  to  23.378  and  19.072“C.  Since  M is  assuaed 
to  be  zero  at  the  ocean  surface,  equation  (10)  can  now  be  solved  for 


-146  gpj/cm'^/hr,  which,  for  the  colixm  one 


square  centimeter  in  cross-section  corresponds,  since  p was  fixed  at  a 
value  of  one  in  equation  (10a),  to  a downward  velocity  of  1.46m/hr.  This 
is  the  downward  velocity  necessary  to  advect  the  heat  to  produce  the  change 
observed  in  the  logarithmic  tenperature  profiles.  For  the  mass  transport 
values  at  the  200m  and  all  deeper  levels,  the  value  of  in  equation 
(10)  is  zero.  The  value  for  M in  the  100 — 200m  layer  is  sinply  the 
valvje  of  M in  t±ie  0 — 100m  l^er  inmediately  above  and  the  T^  is 

obtained  from  equation  (13) . The  ri^t-hand  set  of  arrows  in  Figure  4 
shows  the  values  of  vertical  velocity  calculated  using  equation  (10a)  at 
each  level  for  the  observation  time  mder  oonsideration. 


In  sxnmary,  the  vertical  velocities  at  levels  of  100m  increments 
were  calculated  for  1500  GMT  23  Septenher  1975  two  ways,  first  using 
equation  (5a)  and  then  using  equation  (10a) . 

Only  the  rig^t-hand  set  of  values  in  Figure  4 satisfies  the  heat 
budget,  and  may  be  compared  with  the  left-hand  set  to  estimate  the 
effect  of  omitting  the  horizontal  heat  flux  as  in  equation  (5a) . A 
similar  set  of  calculations,  using  both  die  equation  (5a)  method  and  the 
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equation  (10a)  nethod  was  made  for  each,  of  die  55  observation  times  from 
1300  (JiT,  20  September  to  1200  GMT,  26  Septenber  1975,  vdiich  is  when  the 
meteorological  sensors  on  the  buoy  were  turned  off.  Due  to  the  finite 
differencing  scheme  in  time,  this  period  was  expanded,  for  averaging 
purposes,  to  1230  GMT,  20  Septenber  to  1330  GMT,  26  Septenber  1975— a 
span  of  145  hr. 


III.  RESULTS 

A.  RESULTS  OF  THE  MIXED  LAYER  DEPTH  STUDY 

Table  1 is  a list  of  the  estimated  values  of  MLD  verstis  time  for 
t±e  period  1200  CUT,  20  Septenber  to  1200  QMT,  26  Septenber  1975.  The 
colum  in  Table  1 headed  DTG  gives  the  day  of  Septenber  by  the  first 
two  digits  and  the  GMT  hour  by  the  last  two  digits.  Figure  5 shcx«s  the 
same  variables  graphically  for  the  period  1200  GMT,  20  Septenber  1975 
to  1200  GMT,  11  October  1975.  The  ordinate  in  Figure  5 is  positioned 
at  the  of  passage  of  Eloise  over  EB-10 — 0300  GMT,  23  Septenber  1975. 
Notice  that  according  to  this  graph,  the  average  depth  of  the  mixed 
l^er  was  rou^ily  33  meters  prior  to  storm  ^>proada  and  deepened  grad- 
ually as  the  hurricane  winds  Increased.  An  upwelling  of  the  thexmocline, 
i.e.  a decreasing  of  mixed  layer  depth  began  at  the  time  of  hurricane 
passage.  At  1500  QfT,  23  Septenber  1975,  or  about  12  hr  after  eye 
passage,  the  MLD  readied  a tTrin<ni«ii  of  21.98m  as  calculated  from  the  re- 
gression tenperature  profile  for  this  observation  time.  Follcwing  that 
time  the  MLD  swung  throu^  at  least  three  distinct  oscillations  with, 
nearly  perfect  inertial  cycle  periodicity  ( the  inertial  period  at  the 
latitude  of  EB-10  is  26.0  hr).  These  oscillations  are  superposed  on  a 
continued  deepening  trend  before  finding  an  equilibrlun  vacillation 
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with  a CGchined  inertial-dlumal  period  centered  around  a depth  of  about 

52m.  Ihxis  the  net  mixed  layer  deepening  due  to  storm  passage  was  about 

19m.  Notice  in  Figure  6 the  sqplltude  of  the  first  three  wave  peaks 

following  passage  of  the  eye  of  the  storm  over  EB>10.  All  three  of 

these  oscillations  of  the  thermocllne  are  slgilflcantly  larger  in  anpll- 

tude  dian  those  either  before  or  after  the  storm.  If  the  average  anpU- 

tude  of  the  post-storm  oscillations  is  sihtracted  from  t±ie  aoplitude  of 

the  three  largest,  the  anpUtude  of  each  successive  oscillation  decreases 

by  rou^y  e”^  of  the  prior  value.  Also,  as  shown  in  Figure  6,  if  an 

approximation  of  the  zero  crossing  Is  drawn  throu^  the  osclUaticns , a 
■*lct 

rou^  curve  of  e is  produced.  The  rate  of  dampening  of  these  large 
oscillations  is  irportant  for  predicting  the  duration  of  the  storm  wake 
and  for  verification  of  nuoerlcal  models. 

B.  BESOLTS  OF  THE  ISOTHERM  POSITION  STUDY 

Figure  7 is  a graph  of  seven  isotherms  within  the  thermocllne  during 
the  three  days  before  and  18  days  after  Hurricane  Elolse  passed.  In 
ccnparlng  this  gr£ph  with  Figure  1 — NDBO's  version  assuning  linear 
tenperature  gradients — It  can  be  seen  that  overall,  the  general  trend 
of  the  ipper  isotherms  is  the  same , each  graph  having  three  distinct 
oscillations  followed  by  a net  deepening  after  hurricane  passage.  One 
notes  In  Figure  7 that  the  Isodierm  depths  are  relatively  uniform  prior 
to  hurricane  passage.  The  first  three  upwelllng  cycles  are  distinct 
and  nearly  uniform  throug^iout  the  ipper  thermocllne.  Following  this 
period  the  oscillations  are  much  smaller  in  the  upper  thermocllne.  In 
agreenent  wldx  the  estimated  mixed  layer  depths  in  Figure  5.  However, 
die  oedllatlons  continue  with  considerable  atqplltude  at  greater  depths. 
It  Is  iaportAt  to  note  that  over  the  entire  period  the  mean  depth  of 
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DTG 

MLD(;m) 

DTG 

MLDCm) 

2012 

31.03 

2309 

36.44 

2013 

33.26 

2312 

11. kl 

2015 

33.68 

2315 

21.98 

2018 

30.12 

2318 

26.39 

2021 

31.74 

2400 

41.90 

2100 

34.14 

2403 

49.31 

2103 

33.27 

2406 

50.14 

2106 

33.59 

2409 

48.90 

2109 

32.95 

2412 

43.83 

2112 

31.70 

2415 

35.05 

2115 

33.75 

2418 

35.41 

2118 

33.05 

2421 

40.42 

2121 

33.41 

2500 

44.91 

2200 

34.27 

2503 

52.02 

2203 

34.20 

2506 

53.67 

2206 

37.55 

2507 

53.74 

2209 

35.60 

2508 

53.56 

2212 

36.60 

2509 

53.16 

2215 

35.47 

2510 

52.30 

2218 

35.29 

2511 

51.20 

2221 

35.77 

2512 

50.12 

2222 

35.35 

2515 

46.20 

j/  • 10 

lO 

2300 

38.47 

2521 

42.24 

2301 

39.31 

2600 

43.37 

2303 

41.85 

2603 

46.92 

2306 

39.08 

2609 

54.04 

2612 

53.34 

tabt-F.  1,  Depth  of  mixed  layer  with  time.  Left-hand  colum  gives  the 
j day  of  Septeober  by  the  first  two  digits  and  the  GMT  hour  by  the  last 

^ two  digits . Values  of  mixed  layer  depth  in  the  ri^t-hand  colvun  were 

fomd  by  entering  the  respective  regression  equation  with  the  surface 

JteD|)erature  to  get  the  presstxre  at  the  bottom  of  the  mixed  layer,  and 
then  converting  from  pressure  to  depth  in  meters. 

.'i 

i 


FIGURE  5.  Graph  of  mixed  la^er  depth  vs.  time.  Ordinate  shows  depth 
in  meters  and  is  positioned  along  x-axls  at  time  of  hurricane  passage. 
The  average  mixed  layer  depdi  was  rou^ly  33ni  prior  to  stoxm  approadi. 
At  hurricane  passage  a 12hr  period  of  upwelling  began,  followed  by 
three  distinct  oscillations  of  the  thermocline.  Average  depth  of  the 
mixed  layer  after  these  three  distinct  oscillations  was  rot^ly  52m. 
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FIGURE  6.  Decay  of  intaxnal  wa^a  aopUcude  and  mixed  layer  deepezdng 
rata.  A roug^  approodoatlGn  of  dia  average  post-storm  aoplltude  of 
diermocline  oscillation  was  s\l>tracted  frm  die  aopUtude  of  die  three 
distinct  oedllatlana  imnedlately  after  hurricane  passage.  The  re- 
maining length  of  the  downward  strokes , indicated  by  the  three  axrom , 
was  th«  measiired.  The  values  were  foend  to  decrease  by  rou^ly  e~^ 
of  die  former  value.  Also  a rou^  approximation  of  the  zero  crossing 
(heavy  blade  line)  sxiggests  diet  the  rate  of  deepening  of  die  mixed 
layer  is  e"*^^. 
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each  isothezm  does  not  change  significantly,  vAildi  Is  equivalent  to 
stating  that  the  mean  temperature  in  the  colunn  does  not  change  signifi- 
cantly. Nevertheless,  there  Is  a net  depression  of  the  Isotherms  for 
£^xxit  10  days  following  hurricane  passage.  In  this  perlcxl  there  must  be 
a net  downward  vertical  motion.  Further  investigation  of  these  Isotherms 
and  the  meaning  of  their  hanging  position  in  time  applies  to  the  prob- 
lem of  balancing  the  heat  budget,  the  resiolts  of  hlh  are  dlscnjssed 
next. 


C.  RESULTS  OF  THE  HEAT  BUDGET  STUDY 

Figure  8 shows  the  calculated  logarithmic  tenperature  profiles  for  a 
pre-storm  (solid  line)  and  a typical  post-storm  (dashed  line)  observa- 
tion Hmo  down  to  the  200m  level.  As  expected,  there  was  a cooling  and 
deepening  of  the  mixed  layer.  The  warming  of  the  water  colucm  below  the 
mixed  layer  supports  the  net  Isotherm  lowering  observable  in  Figure  7 
and  Is  suggestive  of  heat  transport  by  downward  vertical  velocities . It 
was  desired  to  solve  the  heat  budget  to  calculate  these  vertical  velo- 
cities. But  an  oscdllating  thennocline  req[ulres  averaging  vertical 
velocity  over  an  integer  nmber  of  cycles,  otherwise  the  heat  budget  Is 
obscured  by  apparent  heat  storage.  Hhce  It  was  decided  to  solve  the 
heat  budget  in  four  parts  (see  Figure  9) , the  first  part  being  the 
nearly  three-day  period  of  storm  approach,  the  second  part  being  the 
perlcxl  of  tremendous  ipuelllng  from  the  time  of  hurricane  passage  xzitll 
1330  OCT,  23  Septaober,  the  third  part  being  the  period  fron  13X  GMT, 

23  Septcober  to  16X  GMT,  25  Septanber  (a  51hr  period  at  a latitude 
where  inertial  cycles  arc  26.0hr,  chosen  because  of  the  distinctiveness 
of  the  periodicity) , and  die  fourdi  part  being  the  period  from  16X  OfT, 
25  Septaober  until  13X  (KT,  26  Septaober— one  and  one-half  hours  after 
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FIGURE  8.  Graph  of  teaperature  vs.  depth  for  typical  pre-storm  (solid  line)  and  post-storm 
(dashed  line)  cbservaticn  times.  MLD  deepened  fi'cm  rou^ly  33m  to  52m  vhile  cooling  more 
than  2®C.  Notice  isotherm  lowering. 


Time  (hours) 


th«  meteorological  sensors  on  the  b\x}y  were  tiimed  off.  The  total  length 
of  averaging  was  145hr.  Vertical  velocities  for  each  of  the  55  ob- 
servation times  during  this  145hr  period  were  calculated  two  ways,  using 
equation  (5a) , and  using  equation  (10a) . Table  2 shows  the  valxjes  ob- 
tained through  equation  (5a)  under  the  colums  headed  w and  the  values 
obtained  through  equation  (10a)  under  the  colums  headed  . The 

values  are  given  for  five  depths  for  all  55  observation  times  considered. 

The  mean  values  for  each  of  the  four  parts  of  the  heat  budget  were  cal- 
culated by  wel^ting  each  value  with  the  appropriate  time  interval  in 
the  far- right  colum,  which  is  one-half  of  the  time  span  between  the 
previous  and  the  subsequent  observations. 

Figure  10  shows  schematically  the  mean  values  over  the  entire  peri.od 
listed  in  Table  2 for  the  five  levels  considered.  The  corresponding 
horizontal  mass  transports  were  obtained  using  ecpjatlcn  (§) . Individual 
values  of  w vs.  at  all  depths  in  Table  2 together  with  the  in- 

formation supplied  by  Figure  10  Illustrate  the  following  points  regarding 
the  four  parts  of  the  heat  budget. 

; 1.  Part  One 

t During  this  period  the  vertical  velocities  deduced  by  both 

i methods  are  generally  similar  in  sign  and  magnitude.  Larger  vertical 

motions  are  found  by  both  methods  during  the  24  hours  prior  to  the  end 
of  the  period  (time  of  hurricane  passage) , particularly  at  lorer  levels . 

t 

These  values  are  several  times  the  mean  values  over  the  period,  vhldi  | 

s 

are  generally  downward.  | 

2.  Part  TWO 

This  period  was  selected  to  Illustrate  the  ri^ld  upwelling  Just 
after  hurricane  passage.  Average  values  of  Sm/hr  for  this  12 -hour 
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period  are  not  tncccmon.  The  w values  frcm  equatlai  (5a)  without 
horlzcxital  teqperature  advectlon  effects  Incliaded  have  a more  uniform 
vertical  profile  (left  side  of  Figure  10,  Part  2)  than  those  calculated 
frcm  equation  (10a).  However,  as  noted  before,  the  values  on  the  left 
side  of  Figure  10  do  not  satisfy  the  heat  budget. 

3.  Part  Three 

During  this  period  the  thermocline  region  is  oscillating  from 
its  hipest  u{ward  excursion.  It  may  be  noted  in  Table  2 Part  3 that 
the  positive  and  negative  vertical  notions  are  not  in  phase  after  about 
0300  CUT,  24  Septenber  1975.  As  indicated  in  Figure  10  Part  3,  the 
lower  level  mean  vertical  motion  becomes  \^3ward  at  500m  whereas  the 
levels  nearer  the  surface  are  still  experiencing  downward  vertical  motion. 

4.  Part  Four 

This  portion  of  the  record  contains  the  third  oscillation  of  the 
thermocline  following  the  hurricane  passage.  It  appears  that  the  ver- 
tical motion  oscillation  at  100m  is  about  90  degrees  out  of  phase  with 
that  at  500m  by  0900  Off,  26  Septenber  1975.  Rather  large  mean  verti- 
cal motlan  values  (Figure  10  Part  4)  are  shown  at  several  levels  be- 
cause not  a cooplete  inertial  cycle  is  inclxxled  in  the  period.  The 
vertical  structure  of  alternating  inflow  and  outflow  with  depth  on  the 
ri^t  side  of  Figure  10  Part  4 contrasts  with  the  more  uniformly 
varying  Inflow/outflow  vertical  structure  deduced  from  vertical  motion 
as  calculated  via  equation  (5a) . 

D.  EBBDR  DISCUSSION 

According  to  Uithee  and  Johnson  [1976] , the  data  were  accurate  to 
2 

within  O.lk^Vcm  in  presstire,  representing  ^out  one  meter  in  depth,  and 
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FIGURE  10.  Mean  values  of  vertical  velocities  at  five  depths  (four- 
part  heat  budget) . Part  1 is  on  this  page  and  Parts  2 , 3 and  4 are 
on  the  next  three  pages.  Valxies  of  w in  left  colxnn  were  obtained 
using  equation  (Sa),  i.e.  assuning  zero  horizontal  teoperature  advec- 
ticn.  Values  of  M /p  in  ri^t  coliim  were  obtained  using  equation 
(10a).  i.e.  not  assdinhig  zero  horizontal  temperature  advection.  The 
corresponding  values  of  horizontal  mass  traisport  were  obtained  using 
equation  (8) , represent  the  total  mass  transp^  across  all  vertical 
faces  of  die  colixm.  and  have  the  units  of  hectograms  per  hour. 


iO-OS^C  in  ten|)erat:ure  except  near  the  top  of  the  theimocline  where 
the  accuracy  dropped  to  ±0.20‘*C.  Ideally  then,  the  errors  induced  by 
die  regresslcn  tazperature  profiles  should  fall  within  these  tolerances. 
Tjhle  3 shows  the  observed  and  regression  teoperatures  and  the  error 
for  each  of  the  times  used  in  Parts  2,3  and  4 of  the  heat  budget  calcu- 
lation, all  of  which  were  more  erroneous  than  part  1.  Since  the  colums 
headed  , T2QQ  and  are  the  tesqperatures  observed  by  the  three 

substirface  sensors  on  EB-10,  it  is  emphasixed  that  these  measiirements 
were  as  the  sensors  migrated  above  and  below  these  nominal  depths 
in  response  to  buoy  movement. 

The  root  mean  square  regression  temperature  error  for  the  three 
colums  are  shown  in  Table  3.  Note  that  at  SOm  the  teqperature  ranges 
from  chout  23*C  to  26.7*C,  the  largest  error  in  fitting  the  profile  is 
^xxit  0.25*C.  Uhereas  the  error  values  are  somewhat  outside  the  ex- 
pected accuracy  of  the  instrunents , it  should  be  noted  that  the  storage 
in  ecpation  (10)  depends  more  on  the  tenperature  gradient  than  on 
the  absolute  value  of  taiperature.  The  smooth  variation  in  time  of  the 
temperature  profiles  derived  by  this  medicxl  were  irportant  for  the  heat 
bxxlget  calculations. 

The  accuracy  of  calculating  the  MLD  using  the  logarithmic  tapera- 
ture  profiles  deserves  cannent.  Of  course  die  accuracy  of  the  MLD  is 
limited  by  the  accuracy  of  eadi  regression  curve.  Given  an  accurate 
regression  curve  and  a surface  tuoperature  accurate  to  within  ±0.20*C, 
die  MD  as  calculated  in  diis  thesis  would  be  knoMi  to  within  t2m. 

The  accuracy  of  the  MLD  is  certainly  governed  by  the  accuracy  of  the 
method,  i.e. , the  very  idea  of  calculating  the  MLD  assuning  an  isother- 
mal mixed  layer  above  a logarithmic  tcaperature  profile.  This  model 
has  not  been  used  in  the  past. 


ETC 

^p500 

“®30 

^p200 

°®200 

^p500 

®*300 

2303 

26.727 

-0.118 

16.462 

0.324 

9.810 

-0.205 

2306 

26.916 

-0.109 

17.159 

0.291 

11.027 

-0.182 

2309 

25.887 

0.033 

16.417 

-0.080 

10.094 

0.047 

2312 

23.334 

0.164 

15.009 

-0.376 

9.547 

0.211 

2315 

22.936 

0.247 

14.612 

-0.543 

9.171 

0.296 

2318 

23.332 

0.193 

14.586 

0.439 

8.860 

-0.245 

2A00 

25.549 

-0.124 

15.788 

0.342 

9.348 

-0.217 

2A03 

26.290 

-0.085 

16.282 

0.225 

9.696 

-0.141 

2A06 

26.181 

-0.023 

16.216 

0.058 

9.600 

-0.035 

2A09 

25.776 

-0.088 

15.268 

0.232 

8.309 

-0.144 

2412 

25.160 

0.017 

15.071 

-0.044 

8.400 

0.026 

2415 

24.143 

0.152 

14.842 

-0.353 

8.705 

0.200 

2418 

24.137 

0.036 

14.936 

-0.089 

8.888 

0.053 

2421 

24.799 

0.057 

15.336 

-0.139 

9.157 

0.082 

2500 

25.294 

-0.008 

15.487 

0.021 

9.072 

-0.013 

2503 

26.240 

-0.084 

15.821 

0.224 

8.968 

-0.139 

2506 

26.647 

-0.078 

16.052 

0.205 

9.048 

-0.127 

2507 

26.736 

-0.059 

15.961 

0.154 

8.850 

-0.095 

2508 

26.731 

-0.054 

15.778 

0.141 

8.568 

-0.086 

2509 

26.618 

-0.078 

15.622 

0.206 

8.426 

-0.128 

2510 

26.483 

-0.101 

15.538 

0.272 

8.323 

-0.170 

2511 

26.363 

-0.112 

15.505 

0.299 

8.276 

-0.187 

2512 

26.155 

-0.119 

15.466 

0.320 

8.343 

-0.201 

2515 

25.555 

-0.098 

15.267 

0.261 

8.431 

-0.163 

2518 

24.672 

-0.039 

15.007 

0.102 

8.682 

-0.063 

2521 

25.295 

-0.097 

15.389 

0.263 

8.920 

-0.165 

2600 

25.450 

-0.044 

15.540 

0.114 

9.065 

-0.070 

2603 

25.943 

-0.070 

15.660 

0.184 

3.390 

-0.114 

2606 

26.647 

-0.038 

15.807 

0.100 

8,753 

-0.061 

2609 

26.986 

-0.061 

15.807 

0.160 

8.526 

-0.098 

2612 

26.828 

-0.137 

15.424 

0.378 

7.983 

-0.240 

RMS 

0.102 

0.255 

0.153 

TABLE  3.  Error  in  regresslcn  temperatures  as  ccqpared  to  observed  tem- 
peratures. Error  is  defined  as  regression  teD|)erature  (T  ) minus 

observed  teoperacure  (T)  . Both  T and  T are  for  sensor  ^ 
depdi  (vhidi  varied) , although  the  ^ nominal  depths  of  50  , 200  and 
SOOm  were  used  in  naming  these  variables . Root-mean-square  of  the  three 
j colums  of  errors  is  given. 
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IV.  SlMl^RY  AND  CONCLUSIONS 

Ihe  buoy  Q-10  measured  the  temperature  at  the  sxirface  and  three  time- 
varying  subsxirface  depths  during  die  passage  of  Hurricane  Eloise.  In 
order  to  study  mixed  layer  deepening,  thermocline  response  and  the  heat 
budget,  temperature  profiles  were  needed  that  represented  the  thermal 
structure  of  the  ocean  from  surface  to  depth  following  storm  passage. 

These  profiles  were  calculated  assuning  an  isothermal  mixed  layer  ^ove 
a thermocline  with  temperature  proportional  to  the  natural  logarithm  of 
hydrostatic  pressure.  The  resulting  graj^  of  mixed  layer  depth 
versus  time  showed  that  prior  to  the  arrival  of  Eloise  at  E3-10,  the 
average  mixed  layer  depth  was  about  33m.  As  the  winds  Increased  due  to 
hurricane  approach,  die  mixed  layer  deepened  steadily  to  about  42m  before 
upwelUng  to  approximately  22m.  The  thermocline  then  underwent  three  dis- 
tinctly large  oscillations  of  inertial  periodicity,  ihlle  the  mixed  layer 
continued  to  deepen.  The  post-storm  average  mixed  layer  depth  was  about 
52  meters. 

Vertical  velocities,  calcnilated  first  by  assuning  zero  horizontal 
tenperature  advection  in  the  material  derivative  ecpiatlcxi  and  second  by 
finding  the  mess  transport  necessary  to  balance  the  heat  budget,  show 
diet  in  die  igiper  500m  of  the  water  colum  dowtsrard  vertical  motion  of 
Iffl/hr  or  less  prevailed  during  storm  appraodi,  followed  by  upward  verti- 
cal velocity  as  great  as  5.3Sm/hr  during  the  12  hours  imnedlately  follow^ 
ing  hurricane  passage.  Next,  during  a 51-hour  period  in  whidi  die  ther- 
mocline underwent  die  first  two  of  three  large  oacillatlcna , downward 
vertical  valocity  cn  the  ordar  of  0.5m/hr  prevallad  in  dia  upper  300m  of 
the  water  colum  with  sll^t  upward  valocity  at  tha  40Qm  aid  SOOm  levela . 
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The  heat  budget  could  not  be  solved  during  the  entire  period  of  the  dilrd 
large  oscillation  of  the  themocllne  because  the  meteorological  sensors 
on  buoy  had  been  shut  off,  but  during  the  first  half  of  die  third 
oscillation,  i.e.,  the  downward  stroke,  calculated  vertical  velocity  was 
generally  downward,  as  great  as  1.27m/hr. 

Aldiou^  die  time-averaged  values  of  vertical  velocities  obtained  by 
neglecting  horizontal  tenperature  advectlon  were  generadly  in  good  agree- 
ment with  those  obtained  considering  horizontal  tenperature  advectlon  in 
all  four  parts  of  the  heat  budget  study,  when  cooparing  the  individual 
values  cooprising  the  averages  of  the  two  methods,  one  sees  that  the  mag- 
nitudes of  the  velocities  differ  significantly.  The  values  obtained  by 
considering  horizontal  advectlon  are  assured  to  be  more  accurate. 
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